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Abstract 

The effective Lagrangian and power counting rules for non-relativistic gauge the- 
ories are derived via a systematic expansion in the large c limit. It is shown that 
the 1/c expansion leads to an effective field theory which incorporates a multi- 
pole expansion. Within this theory there is no need for heuristic arguments to 
determine the scalings of operators. After eliminating c from the lowest order 
Lagrangian the states of the theory become independent of c and the scaling of 
an operator is given simply by its overall coefficient. We show how this power 
counting works in the calculation of the Lamb shift within the effective field 
theory formalism. 



Effective field theories are indispensable tools for studying systems with disparate scales. 
The idea dates back to the Euler-Heisenberg Lagrangian for QED [T|] and has been utilized in 
the context of calculating strong interaction corrections to various processes. Recently, the 
application of effective field theories in heavy quark systems has led to great progress in our 
understanding of weak decays. In particular, heavy quark effective field theory (HQET) has 
been utilized to study hadrons composed of one heavy quark. The use of HQET allows us to 
separate the physics stemming from the two scales which are relevant to heavy-light bound 
states, namely the heavy mass m and the strong interaction scale Aqcd- In a seminal paper |§ 
Caswell and Lepage introduced a similar effective field theory to study non-relativistic bound 
states. However, this theory differs from HQET in several very important ways. The de- 
scription of non-relativistic bound states is complicated by existence of the small parameter 
v ~ a(mv) in the effective theoryf], where v is the relative velocity of the particles of mass 
m which compose the bound state. Furthermore, in HQET the heavy quarks are labeled by 
velocities which are unchanged by bound state effects at leading order. Velocity changing 
weak transitions are accounted for by "integrating in" quarks of varying velocities ||. In 
heavy-heavy systems it is no longer true that the quark[| velocity is fixed. Indeed Coulom- 
bic, velocity altering exchanges are what builds up the Schroedinger kernel. Thus, it is clear 
that the effective Lagrangian for the two systems should be dissimilar. Indeed, in HQET 
an operator scales in 1/M according to its dimension, whereas as in non-relativistic gauge 
theories this is not so unless one modifies the Lagrangian, as will be shown below. 

HQET relies upon an expansion in the heavy quark mass, whereas Non-relativistic Gauge 

Effective Field Theories (NRGT) are expansions in the relative velocity, or equivalently, as 

we shall couch it c, the speed of light. Though it is well known that certain identical 

operators of the same mass dimension may be of different orders in their respective expansion 

(most notably ip^(D 2 /2m)if)), it has not been pointed out that even for fixed dimensions 

the operators of the two theories will in general be different. For instance, if one uses 

dimensionally regulates HQET mixing can only occur between operators with the same 

lr This complication is often said to lead to the problem of having many scales in the theory to, mv, mv 2 . . . 
2 We will refer to the bound state constituents as quarks though they may be electrons as well. Further- 
more, we call the gauge particles gluons to generalize to the Non-Abelian case. 
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scalings in 1/M, this will not be true in non-relativistic field theories unless one modifies the 
Lagrangian, as will be shown in this letter. 

In the effective field theory formalism we write down a low energy Lagrangian which 
reproduces the S matrix elements (in some cases the Greens functions) of the full theory up 
to some chosen order in a double expansion in the coupling a and some other parameter 
which dictates the size of the matrix elements in the low energy theory, the heavy quark 
mass in the case of HQET. The difference between the full theory and the effective theory 
lies in the ultraviolet modes. This difference is accounted for in the low energy effective 
theory by the proper choice of coefficients in the Lagrangian (the "matching" procedure). 
The utility of the effective theory lies in the fact that calculations in the effective theory are 
much simpler now that all the short distance physics has been trivialized. 

A crucial part of the matching procedure is the bookkeeping of the expansion parameters. 
For instance, in HQET where the expansion parameters are a s and 1/m, it is possible to 
place all the dependence on the expansion parameters into overall coefficients of operators in 
the Lagrangian. This is helpful for two reasons: In performing the matching, operators which 
are formally of higher order do not contribute to the coefficients of lower order operators. 
Second, the order at which matrix elements of operators enter should be determined by 
the coefficients of the operators. That is, the states should not depend on the expansion 
parameters in such a way that the power counting is jeopardized. For example, in HQET 
the states are independent of m since the lowest order Lagrangian is independent of m and 
the normalization of the states is chosen to be 

(v' | v) = 2v°(2tt) 3 5 3 (v' -v) (1) 

For the case of NRGT's the expansion parameters are a and v, the relative quark veloc- 
ity. One immediately sees that things will be more difficult in this case since the expansion 
parameter v is dimensionless. This complication led to velocity scaling rules, derived via 
heuristic arguments, which assigned powers of v to fields, operators and derivatives 0. A 
simpler bookkeeping method was presented in []7]] where the authors rescale fields and co- 
ordinates by v in such a way as to make explicit the powers of v in the Lagrangian. Here 
we introduce a slightly different approach which follows simply via an expansion in the now 
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dimensionful parameter 1/c. To the extent that the physical system is truly non-relativistic, 
operators with velocity dimension n are of magnitude ~ (v/c) n , with v a dynamically gener- 
ated scale. This is analogous to HQET where operators of mass dimension n are of magnitude 
(Aqcd / w) n . Moreover the 1/c expansion forces one to modify the Lagrangian. 

Let us consider the large c limit of a non-Abelian gauge theory. In this limit the La- 
grangian density is given by|2|], 

2 C C 4 J C 

+ (iD + 51 U + ■ + —L^DV 

y 2m J 2mc Sm 6 c z 

+ ^#(D • E - E • D)^ + • (D x E - E x D)^ + 0(l/c 3 ), (2) 

where 

D = --gA ; D = V - ^A, (3) 
at c 

and ^ is a non-relativistic 2-spinor describing the heavy quark. In addition we have rescaled 
the fermion field by a factor of y/c. For simplicity we have omitted a 2-spinor describing 
the heavy antiquark. The constants Cp, cd and cs are determined by matching onto the full 
theory. 

The explicit powers of the dimensionful parameter 1/c in @ now makes the power 
counting simple. The lowest order Lagrangian is 

C = ^A a - ~d A?) 2 - i(fitAJ - djAf) 2 + ^ (iD + J,. (4) 

Notice that we have retained the 1/c term in the kinetic energy of the transverse gluons 
in the lowest order Lagrangian. That this is necessary is easily seen by considering the 
Hamiltonian, in which the coefficient of the kinetic energy is c 2 . The eigenstates of this 
lowest order Hamiltonian constitute the states of the effective theory. 

At this point we still have not accomplished what we set out to do, namely, trivialize the 
c dependence of the Lagrangian. Eq. (|), as it stands, will lead to a transverse gluon prop- 
agator with non-trivial c dependence which can, and as we shall see below does, jeopardize 
the power counting in 1/c. To fix this problem we have a choice to either rescale the time 
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or spatial coordinates of the gauge field by c. However, rescaling the time coordinate of the 
gauge field is unacceptable as it will destroy the initial value problem because, the Hamilto- 
nian for the gauge field and fermion fields would then depend on different time coordinates. 
Thus, we make the rescaling 



Ai(y = x/c,t) = y/cAi(x,t), 



(5) 



leaving the Coulomb gauge Lagrangian 



<fWW>4 ? ) - + d>x 



a\2 

o) 



(6) 



We choose to work in the Coulomb gauge for the rest of the paper since it is the most natural 
choice in a non-relativistic theory. Moreover, it allows for the clean separation of powers of 
1/c as is clear from (^). Note that the states of the effective theory are the eigenstates 
of the lowest order Hamiltonian derived from this Lagrangian and, as such, are unconfined 
Coulombic bound states. While these states are independent of c they are not independent 
of g. The Coulomb gluons are leading order and are not treated perturbatively. We will 
return to the issue of confining effects at the end of the paper. 

Now if we follow through with the consequences of the 1/c expansion we will be forced to 
incorporate the multipole expansion. To see this, let us now consider the 1/c corrections in 
(0), concentrating for the moment on the Abelian pieces (the extension to the non-Abelian 
case follows trivially) . The leading 1/c corrected fermionic bilinear Lagrangian is given by 



d 3 x^(t,x) 



ip(t,x). 



(7) 



Expanding in 1/c leads to the multipole expansion 



C 



nip 



mc 



3/2 



^(t,x) (A(t,0) + ^V-A(t,0)A + 



^(t,x) 



+ 



CF 



where 



B 



2 mc 5/2 

d 



^(t,x)a ■B(t,0)ip(t,x) + ■ 



d{yi 



;Mt,y). 



(8) 



(9) 
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The expansion breaks translational as well as gauge invariance which are symmetries that 
are restored at each order in 1/c, and the coefficients in front of each operator are therefore 
fixed. 

Now let us consider the matching procedure for this theory. This will elucidate the 
power counting scheme as well as the problems one runs into if the multipole expansion is 
not performed. Using standard diagrammatic techniques it is easy to show that in the full 
theory the one particle irreducible amputated iV point functions with L loops has an overall 
factor of f\ 



Use of dimensional regularization leads to the functional dependence Y N {qi-qj/ (mc) 2 , fi 2 / (mc) 2 , e). 
On shell these corrections will be both IR as well as UV divergent. The UV divergences are 
taken care of by renormalizing the full theory, while the IR divergences will cancel with those 
in the effective theory, since both theories behave the same in the infrared. To match onto 
the effective theory we then expand this full theory result in a power series in 1/c. The 
coefficients q are then chosen so that the above expansion is reproduced by the effective 
theory, which we now discuss. 

At tree level the matching is trivial. Beyond the terms in (7) the effective theory Hamil- 
tonian will also contain spatially non-local, instantaneous four quark operators which scale 
as 1/c 2 . These operators arise as a consequence of taking the leading term in the c expansion 
of full theory diagrams with transverse gluon exchange between quarks. 

Let us now study the one loop correction in the effective theory. We will first perform 

the calculation using (fj) and show that operators which are supposedly of higher order in 

1/c will generate lower order operators, even within dimensional regularization. We will 

then show that using (0) no such mixing occurs, even within dimensional regularization. 

Let us consider the calculation of the correction to the two point function coming from 

two insertions of the magnetic moment operator, a ■ B. If one wishes to keep the power 

counting such that a given operator should scale as a fixed power in 1/c (as dictated by its 

overall coefficient), the higher order operators will not contribute to the renormalization of 

3 Since a ~ e 2 /c we could equally write the expansion purely in terms of 1/c taking e w 1. We have 
chosen to write the expansion in terms of a in analogy with the expansion in HQET. 




(10) 



5 



the lowest order Lagrangian. That this is so in heavy quark effective theory is easy to see 
simply on dimensional grounds. 

Using the Feynman rules derived from (Q) we have 

(2) _ c 2 F C(R)e 2 r d n k (k x a) ■ (k x a) 

4m 2 c 2 J (2tt)" (k 2 /c 2 -k 2 + ie)(E + k - (p + k) 2 /2m + ie)' 

We integrate over ko and choose to close the contour in the upper half plane, picking up only 
the negative energy pole. Integrating over the magnitude of k leaves 



r (2) = 4C(fl)<* r(2c) 
7T 



2 (E-^)+4mc 2 + -^- 
v 2m' 32m 



+ finite, (12) 



where C(R) is 1 and 4/3 in QED and QCD, respectively. We see that using (|2|) leads to 
the mixing of operators of different orders in 1/c. This can be avoided by expanding (|TT|) in 
powers of k. In previous matching calculations |], |Tl], ^ this is in fact what has been done, 
and it is justified as a bona-fide approximation, as an expansion in small A;||. This amounts 
to dropping the k dependent terms in the denominators. Such an expansion makes working 
with dimensional regularization particularly simple, since we now have a scaleless integral 
which vanishes in this scheme. 

However, we emphasize that, once we choose to amend the Feynman rules we must also 
amend the low energy theory. Not to do so would destroy the power counting scheme. In 
Ref. 0, the authors point out that the transverse gluons can lead to enhancements in c 
(their 1/v ) in the low energy theory. This is only true if one insists upon calculating the low 
energy matrix elements using (Q). Instead one must calculate using an amended Lagrangian 
(8) which reproduces the Feynman rules utilized in the matching calculations. Indeed it is 
simple to show that Feynman rules in (|8D leads to the necessary expansion. 

Using (|SD we may now calculate anew the contribution to the two point function from 
two insertions of a ■ B, 

•p(2) __ c 2 F C(R)e 2 r d n k k? 

2m 2 c 5 J (27i) n (k 2 -k 2 + ie)(E-p 2 /2m + k + ie)' 

Notice that all the c dependence is now explicit, and the variable k has units of energy. A 
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simple calculation leads to the result 



r(2) = 4^ r(2£) (^-£) 3 + flnite (14) 

This correction renormalizes some higher order operator which vanishes by the equations of 
motion. It is clear that once we use the correct effective, theory the corrections resulting 
from the insertion of higher order operators will never feed down into the matching for lower 
order operators. 

Just as the 1/c expansion dictates the proper matching procedure, it also greatly simplifies 
the power counting rules in the low energy theory. In HQET the magnitude of a matrix 
element is dictated by the explicit power of l/m in its coefficient. Thus, an operator of 
mass dimension d with a factor of m~ n in its coefficient is of order QC ° . In NRGT's, as 
formulated here, the power counting is just as simple. All the powers of 1/c have been made 
explicit, and we can read off the size of a matrix element simply by counting powers of c and 
doing dimensional analysis. 

As an example of the power counting procedure let us consider the pedagogical example 
of the relativistic corrections to the bound state energy of Hydrogen] IT]. Given that a and 
v/c are of the same order, we must calculate the matching corrections to the appropriate 
order in a for the accuracy we wish to attain. The leading corrections come from pure (v/c) 2 
corrections stemming from the effective Lagrangian. The first such term is the correction to 
the kinetic energy tp > d 4 /( 8m 3 c 2 )ip. By dimensional analysis, its matrix element is mv 2 (v/c) 2 , 
which yields a relative contribution of order (v/c) 2 . There are no a corrections in the 
matching to this operator. Next there are (v/c) 2 corrections coming from the Darwin term 
and the spin orbit coupling which are again of relative order (v/c) 2 , since we may pick out 
the Coulombic piece of the electric field. Both the Darwin term and the spin orbit term will 
get matching corrections at order a and will thus give a contribution to the energy shift at 
relative order (v/c) 3 as well. There are no further corrections at order (v/c) 2 , assuming the 
proton to be infinitely heavy so that the magnetic interactions become irrelevant. There is 
a correction of relative order (v/c) 3 coming from corrections to the Coulomb potential due 
to pair creation which are accounted for in the matching and lead to a term in the effective 
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Lagrangian given by 

diAo&diAo . 

O v = c v a t—: . (15) 

m z c z 

Finally, we come to the energy shift due to transverse photon propagation in the bound state. 
This is a self-energy correction of the electron propagating in the Coulomb background and 
is of relative size (v/c) 3 due to a factor of c _3//2 coming from each transverse photon vertex. 
This agrees with the well known result for the Lamb shift. 

Let us now return to the issue of confinement in the Non-Abelian theory. While it may 
be surprising that the non-Abelian couplings are subleading, it is clear that the non-Abelian 
nature of the theory should be irrelevant to the details of the bound state as its size is 
reduced. The confining effects in a Coulombic bound state should be suppressed by powers 
of h.QCD/'m- Setting r ~ 1/mv in the virial theorem 

a s (l/r) mv 2 

, (16) 

r c 

gives v/c ~ a s (mv), which leads to the conclusion that, at small v, Aqcd/^i scales likefj v/c. 
The question then becomes how do we properly take into account the effects of confinement 
in this effective field theory? Since confinement will not arise in perturbation theory we 
must insist that the zeroth order states contain the confining potential. The omission of 
the non-perturbative effects will lead to the breakdown of the 1/c expansion as we will now 
show. 

If we include the full non-Abelian gluondynamics in the zeroth order effective theory, then 
the Coulombic potential will now be modified by the linear rise due to confinement. At first 
this may seem a bit worrisome since the non-Abelian piece has explicit factors of 1/c in front 
of them and thus the states of the theory will depend upon c, which is exactly what we were 
trying to avoid. However, this dependence on c will not destroy the systematics. The analytic 
dependence now implicit in the matrix elements will clearly not upset the systematics of the 
1/c expansion. The non-analytic dependence on c will be introduced through factors such 
as 

^l? = (±QOD\ P (17) 
\ m J 



We will take Kqcd to have the units of mass. 
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where P is some positive power. Thus, the c dependence of the matrix elements due to 
non-perturbative effects can only introduce higher order corrections. If we chose to use the 
Coulombic states, we would never see these higher order effects, and since Aqcd/™> ~ v/c, 
we would possibly miss effects of the same order we wish to keep. 

We have shown that Non-Relativistic Gauge Theory (NRGT) is conveniently organized 
in temrs of an expansion in the dimensionful parameter 1/c. Furthermore, if dimensional 
regularization is used, the organization of the expansion is very simple, there is no operator 
mixing across different orders in the expansion and the order of the operators in the expansion 
can be read off directly from their dimensions without resort to heuristic arguments. We 
emphasize, however, that dimensional regularization is not mandatory. Indeed one may 
formulate NRGT's with, for example, a momentum cut-off A. If A m, with m the quark 
mass, then the multipole expansion is automatic. However, one then has a triple expansion, 
in a, 1/c and A/m. It is not necesary to choose A < m; with A ^> m one can work with a 
double expansion, in a and 1/c only.0 In either case operators of different orders in the 1/c 
expansion mix. The organization by 1/c is still useful because once the matching procedure 
has been completed the order at which an operator enters can still be trivially obtained from 
its velocity dimension. 
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5 In fact, it is unnccesary to perform a multipole expansion, if one is willing to fine tune the coefficients 
in the expansion order by order in perturbation theory. 
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